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Abstract 

The theoretical analysis of p—A0p—A correlation function in 10% most central Au+Au collisions 
at RHIC energy ^/snn = 200 GeV shows that the contribution of residual correlations is the 
necessary factor to obtain a satisfactory description of the experimental data. A neglecting of the 
residual correlation effect, leads to unrealistically low source radius, about 2 times smaller than the 
corresponding value for p—A0p—A case, when one fits the experimental correlation function within 
Lednicky-Lyuboshitz analytical model. Recently an approach accounting effectively for residual 
correlations for the baryon-antibaryon correlation function was proposed, and a good RHIC data 
description was reached with the source radius extracted from the hydrokinetic model (HKM). 
The p — A scattering length, as well as the parameters characterizing the residual correlation effect 
— annihilation dip amplitude and its inverse width — were extracted from the corresponding ht. 
In this paper we use these extracted values and simulated in HKM source functions for Pb0Pb 
collisions at the LHC energy U^nn = 2.76 TeV to predict the corresponding pA and pA correlation 
functions. 
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INTRODUCTION 


The study of pA correlation functions (CF) (along with the study of other two-particle 
correlations such as tttt, KK) allows one to obtain the information about the character 
of evolution of the matter formed in relativistic nuclear collisions, in particular about the 
character of collective flow. 

Another great possibility which is open for a researcher in baryon-(anti)baryon correlation 
analysis is the study of strong interaction between the particles of different sorts with Final 
State Interaction (FSI) correlation technique [l-3|. Since LHC produces copiously hadrons 
of different species, including multi-strange, charmed and beauty ones, the advantage of 
this approach is the ability to analyze the interactions even in exotic particle pairs, hardly 
achieved by other means. 

Fitting the experimental correlation function with some analytical formula, e. g. Led- 
nicky-Lyuboshitz model [l|, allows one to extract the quantities describing both the inter¬ 
action in corresponding pairs and the particle emission region size. Of course, if the source 
function ^(r*), describing the spatial structure of the pair emission, is known, an extrac¬ 
tion of unknown interaction parameters becomes more reliable. The corresponding spatial 
structure can be obtained from realistic collision model that simulates the evolution of the 
system formed in high energy heavy ion collisions. In the recent paper jj| the hybrid variant 
of the hydrokinetic model (HKM) jh-T] was used for this purpose. The choice of HKM is 
highly reliable since it is known to provide a successful simultaneous description of a wide 
class of bulk observables in nuclear collision experiments at RHIC and LHC {q]. The model 
also reproduces well 0 the source functions for pion and kaon pairs in Au-|-Au collisions 
at the top RHIC energy [n|, including non-Gaussian tails observed in certain experimental 
source function projections. 

In in order to extract unknown pA scattering length, the corresponding experimental 
correlation function for 10% most central Au-|-Au collisions at top RHIC energy, measured by 
STAR Collaboration js], was htted with the Lednicky-Lyuboshitz formula using the effective 
source radius tq extracted from the HKM source function. It was also found that ro values 
obtained in HKM for baryon-baryon and baryon-antibaryon correlations are expectedly close. 
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while in the STAR experimental analysis js], where the source radii were considered as free 
£t parameters, the extracted pA radius was ~ 2 times smaller than the pA one. One can 
assume that this apparent difference is due to neglect of residual correlations in js]. Such 
correlations can exist between secondary protons and lambdas if their parents were correlated 
(or the parent of one particle was correlated with another particle in the pair). For taking 
into account the residual correlation contribution to the baryon-antibaryon CF, a modihed 
analytical formula was introduced in jd]. As a result, a good description of the experimental 
pA correlation function is obtained and the pA spin-averaged scattering length is extracted 
from the corresponding £t. 

In the present paper we are going to apply the method developed in for description 
of the RHIC data to predict the pA and pA correlation functions in the 5% most central 
Pb-|-Pb collisions at the LHC energy = 2.76 TeV, utilizing the results of Ref. jd] as a 

starting point. 


II. FORMALISM 


a 


The experimental correlation functions presented in j8| are purity corrected. They are 
obtained from the measured ones as 


(L*) _ ^meas(fe ) 1 . 
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where Cmeas(k*) and Ccorr(k*) are the measured and the corrected CF respectively, and A(fc*) 
is pair purity. The latter is defined as the fraction of pairs consisting of primary, correctly 
identified particles. 


with Lednicky-Lyuboshitz analytical formula 
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where Fi{z) = /z and F 2 {z) = (1 — e~^^)lz. This formula is derived starting 

from the basic equation C{k*) = ^|\l/;^j^*(r*)|^^, where the wave function represents 
the stationary solution of the scattering problem with the opposite sign of the vector k*. 
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The angle brackets mean averaging over the total spin S and the distribntion of the relative 
distances «S'(r*). Since typically the source radius can be considered much larger than the 
range of the strong interaction potential, ^le approximated at small k* by the 

s-wave solution in the outer region; 
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The scattering amplitude f^{k*) here is taken in the effective range approximation 

/y^*) = {-js + - *•) . ('‘) 

where is the scattering length and dg is the effective radius for a given total spin S' = 1 or 
S' = 0. The singlet and triplet weights pi for unpolarized particles (supposing the polarization 
P = 0) are in po = 1/4(1 — P^) = 1/4 and pi = 1/^3 + P^) = 3/4 correspondingly. 

As for the baryon-antibaryon case, following j4|, we fit the experimentally measured 
correlation function Cuncorr{k*) with the following analytical expression 


CuncoAk*) = 1 + A(r)(C(r) - 1) + a{k*){Cres{k*) - 1), (5) 


where X{k*) is purity or the fraction of correctly identified pairs consisting of primary par¬ 
ticles, C{k*) is “true” correlation function approximated by Eq. ([2]), a{k*) is the fraction of 
secondary particles which are residually correlated, a{k*) = a{l — X{k*)) and CrAk*) is the 


residual correlation contribution. The latter is taken in the Gaussian form 


m 


Cres{k*) = 1 - /3e 




( 6 ) 


where /? = A > 0 is the annihilation (wide) dip amplitude and P <C rg is the dip inverse 
width. Since a and (5 enter (I5l) only as a product d/3, the latter is treated as a single 
parameter f3 at fitting. 

The source radii tq in both pA and pA cases are extracted from the Gaussian ht SfnA) = 

_,*2 

(2-y7rro)“^e "‘''o to the angle averaged source function. 
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S{r*) = l/(47r) 
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calculated in the hydrokinetic model js-?]. The latter simulates the evolution of the matter 
in relativistic nuclear collision as consisting of two stages — the hydrodynamic expansion of 
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the matter being in local thermal and chemical equilibrium and gradual system decoupling, 
beginning when the equilibrium is lost. The first stage is described within ideal hydrody¬ 
namics and for the second one the hydrokinetic approach is utilized, based on the Boltzmann 
equations in the integral form, with switching to UrQMD cascade at a space-like hypersur¬ 
face. In current study HKM is taken in its simplified hybrid form {q] with sudden switch 
from hydrodynamic evolution to the cascade at the hadronization hypersurface dehned by 
the isotherm T = 165 MeV. 

The model output consists of generated particle momenta and coordinates, that are 
further used to build different observables. The considered angle-averaged source function 
histograms are hlled using the following procedure (here r* is the particle spatial separation 
in the pair rest frame) 


Y.n=i + r4)/(47r(r*. - r*„jA)] 
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Here r*n and r*™ are the pair rest frame r-coordinates of particles 1 and 2 produced in the 
n—th event, is the r-coordinate of the k-th histogram bin center, the function Sa(x) = 1 
if |a;| < A/2 and 0 otherwise, and A is the size of the histogram bin. 


III. RESULTS AND DISCUSSION 


The initial conditions (IC) for HKM calculations simulating the considered case of 5% 
most central Pb-|-Pb collisions at the LHC energy sJsnn = 2.76 TeV are described in detail 


in jg. 


14| . We assume longitudinal boost invariance, so that the IC are specihed in the 


plane transverse to the beam axis only. The transverse energy density profile ej(r'r) at the 
starting time r, = 0.1 fm/c corresponds to Monte Carlo Glauber model and is calculated in 
GLISSANDO code 15|, where the overall scale factor eo, being the maximal initial energy 
density, is fixed basing on the experimental mean charged particle multiplicity. The initial 
transverse flow in present calculations is absent. 

In Fig. [Done can see the angle averaged pA source function S{r*) obtained in hydrokinetic 
model for considered LHC collisions together with the Gaussian fit to it. The source radius 
value extracted from this ht is tq = 3.76 fm, which is about 1.15 times larger than for the 
RHIC case. The pA source function fitting results in the same source radius value. 
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FIG. 1. The pK angle averaged source function calculated in HKM (markers) and the Gaussian 
fit to it (line). The calculations correspond to 5% most central Pb+Pb collisions at LHG energy 
■s/snn = 2.76 TeV, in pseudorapidity range \t]\ < 0.8. Proton pT range is 0.7 < pT < 4 GeV/c, 
and for lambdas 0.7 < pT < 5 GeV/c. The extracted vq = 3.76 fm. 


Having obtained vq from HKM and fixing /jf and dg according to the paper 13| for 
baryon-baryon case, we can model the corresponding correlation function (see Fig. [2]). As 
compared to RHIC, the LHC pA CF is slightly more narrow (the two corresponding Gaussian 
widths differ by a factor of ~ 1.13) and has lower intercept of about 1.7. 

To build the corresponding baryon-antibaryon correlation function one should determine 
the values of parameters entering ([5]) and ([6]) for the LHC case. As well as for RHIC we 
assume /^ = /* = / and dg = dg = 0. The source radius is again fixed from HKM 
calculation, rg = 3.76 fm. As for the real and imaginary parts of the scattering length, 
Re/g and Im/g, they characterize antiproton-lambda strong interaction and hence are not 
changed when switching from RHIC to LHC. So we can use the values extracted from the 
ht to RHIC pA correlation function shown in Fig. 6 of the paper jd], where the Gaussian 
parametrization ([6]) for the residual correlation contribution C^esik*) is applied, Re/o = 
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FIG. 2. The HKM prediction for purity corrected p — A (B p — A correlation function in the 
LHC Pb+Pb collisions at y'SNN = 2.76 TeV, c = 0 — 5%, \r]\ < 0.8, with 0.7 < pT < 4 GeV/c for 
protons and 0.7 < px < 5 GeV/c for lambdas (red line). The LHC source radius value calculated in 
HKM is ro = 3.76 fm. The Lednicky-Lyuboshitz fit to the top RHIC energy correlation function, 
corresponding to the STAR experiment with rg = 3.23 fm extracted from the HKM source 
function is presented for comparison (blue line). 

0.14 ± 0.66 fm and Im/o = 1.53 ± 1.31 fm. The parameter (3 = a[3 describes the strength 
of residual correlations and the fraction of residually correlated non-primary particles. So, 
this parameter depends on particle interaction kinematics. Thus, it also can be taken the 
same as for RHIC, (3 = 0.034 ± 0.005. As for the R parameter, being some effective size 
associated with residual correlations (Rrhic = 0.48 ± 0.05 fm), one can suppose that for 
LHC it will be larger than for RHIC, approximately proportionally to the source radii 
ratio 1.15, i. e. Rlhc = Such an assumption gives 

-Rlhc = 0.55 ± 0.06 fniQ. Note, that since the mentioned parameter values have errors, the 

^ The calculations show that the model correlation function depends weakly on R changing by such a close 
to unity factor. 
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LHC correlation function cannot be predicted exactly, but only up to uncertainties caused 
by these errors in the £t parameters. The uncertainty in the predicted correlation function is 
calculated as AC„ 


= \ Eh 


dCu 




dxi J ^ dxi dxj 

where Xi are the 4 £t parameters, Re/o, Im/o, (3 and R, ai = Axi are the corresponding 

parameter standard errors, and aij are the covariances of parameters Xi and Xj. 

The determination of purity X{k*) at the LHC is not so explicit, since it depends not 
only on the fraction of secondary pairs, coming from resonance decays, but also on the 
experimental setup, or more precisely, on the fraction of misidentified pairs. To clarify this 
issue we have calculated the fractions of pairs made by particles having different origination 
(primary or coming from certain decays) in the hydrokinetic model for RHIC and LHC 
collisions. The results obtained in both cases are quite close and are presented in Table 1. 
Comparing these fractions with those in Table HI for RHIC Sj, one can see that experimental 
fraction of Pprim — Ap^im pairs A = 0.15 is about 2.5 times lower than in HKM, likely due 
to misidentihcation problem, which takes place in the experiment. In HKM simulations, 
on the contrary, all the produced particles are correctly identihed, that leads to such a 
difference between corresponding fraction values. However, since HKM is a realistic model 
that describes both RHIC and LHC bulk observables well, basing on its results one can 
conclude, that true purities, understood as the fractions of primary pairs, at RHIC and 
LHC should be quite similar. 

In such a situation we demonstrate three different plots for the LHC p — A 0 p — A 
correlation function (see Fig. [SHH]). The hrst one demonstrates the model CF with ALHc(fc*) 
the same as for RHIC 0, ALHc(/i^*) = Arhic(^*)- This function is again more narrow 
than for RHIC. In Fig. |4]one can see our prediction for the LHC correlation function with 
Alhc(^*) = 2.5Arhic(^*)) where the factor 2.5 corresponds to the ratio of corresponding 
primary pairs’ fractions in HKM and in the STAR experiment. As for Fig. [5l it shows the 
purity and residual correlation corrected CFs for LHC and RHIC. They are expressed (as it 
follows from Eq. (|5])) through the uncorrected ones as C{k*) = 1 + (Cuncorr(/i^*) — 3)/\{k*) — 
a{k*){C-[es{k*) — l)/\{k*). The latter “true” function can be easily compared with the 
experimental result (corrected in the same way for purity and residual correlations), since 
it does not depend on the fraction of misidentihed particles in the concrete experiment. As 
compared to RHIC, the LHC curve apparently has smaller amplitude and width. 
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TABLE I. The fractions of pA pairs, primary and coming from different decays calculated in HKM. 
These fractions are quite similar for RHIC and LHC cases. 

IV. CONCLUSIONS 

The first predictions for pA and pA correlation functions, calculated within Lednicky- 
Lyuboshitz and hydrokinetic (HKM) models and accounting for the residual correlation 
effect, are presented for the 5% most central Pb+Pb LHC collisions at the energy y/SNN = 
2.76 TeV. The functions’ behavior is predicted based on the results previously obtained for 
Au+Au collisions at the top RHIC energy. 

Both pA and pA correlation function curves in the LHC case are slightly more narrow 
than those at RHIC. 

The LHC source radii ro, calculated in HKM for baryon-baryon and baryon-antibaryon 
cases are similar, ro = 3.76 fm. They are about 1.15 times larger than the corresponding 
RHIC radii. 

The pair purities (the fractions of pairs consisting of primary particles), calculated in 
HKM for RHIC and LHC are very close and are about 2.5 times larger than the experimental 
primary pairs fraction in the STAR Collaboration experiment at the RHIC. This difference is 
most probably because of particle misidentihcation problem existing in the experiment. That 
is why we present results for purity uncorrected pA correlation functions in two variants: 
with purity A(/c*) similar to that in the experiment at RHIC and with RHIC purity scaled 
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FIG. 3. The same as in Fig. [2] for purity uncorrected p — A (B p — A correlation function. The 
HKM source radius for LHC is rg = 3.76 fm. The purity X{k*) is the same as for RHIC case 8|]. 
The scattering length real and imaginary parts, Re/o and Im/g, are taken from the fit to RHIC 
CF that corresponds to Fig. 6 from where HKM source radius rg = 3.28 fm and the Gaussian 
parametrization Q for the residual correlation contribution C^esik*) are utilized. For the LHC fit 
the Cresik*) parameter (3 coincides with that for RHIC, while parameter R is scaled by the factor 
^LHC/^RHic_ LHC fit is determined up to errors in parameters Re/g, Im/g, (3 and R, that is 
illustrated by the band around the LHC curve. 


by a factor of 2.5, that corresponds to HKM purity. The measured LHC correlation function 
should he somewhere between these two limiting possibilities. Depending on the concrete 
experimental setup and the related fraction of misidentihed particles, the real curve can be 
closer either to one or to another variant. 

We also demonstrate the purity and residual correlation corrected baryon-antibaryon CF 
which should not depend on the experiment details and thus can be easily compared with 
experimental result. This function for LHC has smaller amplitude than for RHIC. 
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FIG. 4. The same as in Fig. [3l but LHC purity X{k*) is scaled as compared to RHIC according to 
primary pA pairs fraction in HKM simulations, where there is no misidentification problem. 
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